Comments on "Single meson production in photon-photon collisions and infrared 

renormalons" |arXiv:0911.5226vl] [hep-ph] 

S. S. Agae\Q 

Institute for Physical Problems, Baku State University 
Z. Khalilov St. 23, Az-1148 Baku, Azerbaijan 
(Dated: February 8, 2010) 

This paper contains our comments on the work of the authors A. I. Ahmadov et al. " Single meson 
production in photon-photon collisions and infrared renormalons" arXiv:0911.5226vl [hep-ph] 27 
Nov 2009]. We draw attention to errors made in this work in calculation of the subprocess cross 
sections within the frozen and running coupling approaches, in deriving the Borel transforms and 
Borel resummed expressions for the cross sections in the running coupling approach. We show that 
>— ' ■ the same fatal errors were made by these authors also in their work published in Phys. Rev. D80, 

016003 (2009). 

<n : 

C\ ■ In the work "Single meson production in photon-photon collisions and infrared renormalons" arXiv:0911.5226vl 
[hep-ph] 27 Nov 2009] [l[ authors A. I. Ahmadov, Coskun Aydin, E. A. Dadashov and Sh. M. Nagiyev in accordance 
Ph ■ with the title of the paper intended to compute the cross section of the single meson M production in the photon- 
qq ' photon collision 7 + 7 — > M + X. 

In the framework of the perturbative QCD (pQCD) this process proceeds through leading and higher twist mech- 
anisms. At the leading twist (LT) order the meson M is produced due to fragmentation of the final quark or gluon 
(appearing in the hard-scattering subprocess) into the meson M. The relevant LT cross section can be computed 
using the standard pQCD methods, applicable to inclusive hadronic processes. The higher twist (HT) mechanism 
Qh| implies production of the meson M directly at the hard-scattering subprocess @. In the HT subprocess emergence 
of the meson in the final state becomes possible in the result of the hard gluon exchange between meson constituents. 
The cross section of the meson production via HT mechanism can be found employing the methods of the pQCD 
elaborated to compute hard exclusive processes and relevant factorization theorems [3H5J . By this way HT corrections 
i~ 1 ' to some inclusive processes were obtained . 

In order to find the amplitude of the HT hard-scattering subprocess one should use the distribution amplitude 
of the meson and perform integrations over the longitudinal momentum fractions x\ and x 2 carried by quark and 
antiquark of the meson. In the case under discussion it takes the following form 
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] where /i 2 R and ^ip are the renormalization and factorization scales, Th(x\, x 2 , fJ-p) - hard-scattering function, 
(%i,%2, ftp) is the distribution amplitude of the meson. One of the important questions to be solved in order 
.{^ to calculate Eq.([T]) is a proper choice for the scales ^x 2 R and y? F . It has been advocated @, E3 that to reduce the 
higher order corrections to a physical quantity and improve the convergence of the corresponding perturbation series, 
the renormalization scale, i.e., the argument of the QCD coupling in a Feynman diagram should be set equal to the 
hard gluon's squared four-momentum. This idea was proposed more than 25 years ago, and used in numerous works. 
Nevertheless, A.I. Ahmadov et al. write: "According to Ref.[26] [A. I. Ahmadov et al., Int. J. Mod. Phys. E15, 
1209 (2006)] should be noted that in PQCD calculations the argument of the running coupling constant in both, the 
renormalization and factorization scale Q 2 should be taken equal to the square of the momentum transfer of a hard 
gluon in a corresponding Feynman diagram" ! 

In exclusive processes and in HT Feynman diagrams describing corrections to inclusive ones, the renormalization 
scale chosen this way inevitably depends on the longitudinal momentum fractions X\ and x 2 carried by the meson's 
quark and/or antiquark fi R ~ xi(x 2 ). Then the pQCD factorization formula (in our case Eq.([T])) diverges, since 
a s ((J, R ) suffers from an end-point 0:1(^2) — >• 0, 1 singularities. Therefore, the HT amplitude may be computed using 
two options: one of them is the standard "frozen" coupling approximation, when one fixes x\ and x 2 equating them 
to their mean values x\ — x 2 — 1/2, and remove a s (fJ,%) as a constant from the integral Eq.([T]). In the second 
approach one allows X\ and x 2 to run in the argument of the a s (fi R ) and calculate the HT amplitude using the 
running coupling (RC) method, and removes divergences appearing in the perturbative expression with the help of 
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the infrared (IR) renormalon calculu s II ill (the Borel transformation, resummed expression for a physical quantity) 
and the principal value prescription [12]. It turns out that this method allows one to estimate power corrections 
arising from the end-point integration regions. 

The RC method was suggested in our works [HI, [TJ] (see also, in order to compute effects of the infrared 

renormalons on the pion and kaon electromagnetic form factors (FFs) . Later the RC method was applied for calculation 
of power suppressed corrections to some exclusive processes, vertex function and HT corrections to semi-inclusive 
processes. Namely, it was used for com put ations of the pion, kaon electromagnetic FFs [HHlllil, 7r°7 [l?} and 
flli V 1 ! electromagnetic transition FFs [18| . for evaluation of the power corrections to the gluon-gluon-7/ vertex 
function [l9[ and to the two pion production in the process 77* — > t ttt f20ll . It was also employed to estimate HT 
correction to the semi-inclusive meson production 7 + h — > M + X [2l| . It is worth noting that the RC method 
deals with power-suppressed corrections to exclusive processes coming from the end-point integration regions. There 
are another source of the infrared renormalon effects in inclusive and exclusive processes, namely one appearing due 
to resummation of diagrams with quark vacuum insertions ("bubble chains") into a gluon line. It is known that 
resummation of infinite number of such diagrams corresponds to computation of one-loop diagrams with the runnning 
coupling a s (— k 2 ) at vertices. The HT ambiguities produced during these calculations were employed to model HT 
corrections to numerous inclusive and exclusive processes (for review, see [ll| ) . In the exclusive processes the coupling 
constant a s runs not only due to loop integration, but also because of the integration in the process amplitude over 
the longitudinal momentum fractions of hadron constituents. Thus the exclusive processes have two independent 
sources of the infrared renormalon effects. 

In Ref. Il| both the frozen and running coupling approaches were used. Employing the RC method elaborated in 
Refs. (T3, Il4 . Il7l - [2"l| . authors of [l[ nevertheless "forgot" to cite them - the real sources used in their calculations and 
in preparation of the text of the paper. Now we want to concentrate on "calculations" performed in Ref. jlj. The 
higher twist subprocesses considered in Ref. jj are 
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In the "frozen" coupling approximation they were calculated by Bagger and Gunion | 
results for the jqi — > Mq 2 subprocess cross section in Eqs. (16) and (17) of Ref. Q: 



The authors presented their 
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In Eqs.([3|) and (|4|) s, u, t are the Mandelstam invariants of the subprocess, e\ and e 2 are the charges of the quark q\ 
and antiquark q 2 . The function Im in Ref. @ is defined as 
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where ,7 Q" V =1 t g 1 and "t g is the average squared momentum transfer carried by the hard gluon in a given subprocess" 
equal to s/2 and —u/2 for two type of relevant Feynman diagrams. 

In our work "Single meson photoproduction and IR renormalons" (2lj | it was demonstrated that namely the HT 
subprocesses contribute to the photoproduction of the single meson 7 + h — > M + X via the HT mechanism. 
Within the RC method the differential cross sections of these HT subprocesses were found in our work 2l|. The 



differential cross section of the subprocess 7^ — > Mq 2l when M is a pseudoscalar or longitudinally polarized vector 
meson, is given by the following expression [Eq.(9) in Ref. |21|]: 
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For the transversely polarized vector meson we get [Eq.(lO) in Ref. [2lj]: 
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In Eqs.([H]),([7]) the functions Iu%\ and -K"i( 2 ) encode an information on the meson distribution amplitude and have the 
forms: 
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Mhi = x 2 s, fj, 2 F1 = s/2; [i R2 = -xi_u, \x 2 F2 = -u/2, 

s, u, t are the Mandelstam invariants of the subprocess, e\ and e 2 are the charges of the quarks q\ and q 2 , respectively. 

Equations ([6]), (O are general expressions valid for mesons with both the symmetric (under exchange xi ■<-> x 2 ) and 
non-symmetric distribution amplitudes. In the frozen coupling approximation, i.e., when 

Mhi = s/2; = -«/ 2 > 

from Eqs. (|6|), (JTJ) for mesons with symmetric distributions the Bagger-Gunion results can be obtained. Indeed, in 
this approximation for the functions /i( 2 ) and ^i( 2 ) we find 
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Here the superscript "0" indicates that the functions are found in the frozen coupling approximation. Having per- 
formed some analytical calculations we get 



dcr H1 (ei,e 2 ) %^ 2 C F a E 1 



9s 2 



-*) 



[eiM/ 2 °-e 2 s^ 2 ]' 



1 1 



(11) 



d(T WJ (ei,e 2 ) _ 647r 2 Ci?a£; -t 



9s 4 



[eiziJ^ - e 2 sK°] 



(12) 



These expressions coincide with the Bagger-Gunion results, if to take into account that in Ref. Q e.\ and e 2 are the 
charges of the quark q\ and antiquark q 2 , whereas in our work they denote the charges of the quarks qi and q 2 . 

In Ref. [l| authors presented results of their calculations of the subprocess jqi — > Mq 2 cross section. In accordance 
with Eqs. (2. 10) and (2.11) of Ref. [l[ these cross sections are given by the following expressions: 
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where 
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and Q\ = s/2, Q\ = -u/2. 

Now using Eq.dH]) [Eq.(2.11) in Ref.Q] for the function D(s,u) we find 
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It is evident that Eqs. ([TUl) and (fT5)) are the subprocess cross sections computed in the framework of the frozen coupling 
approximation. Having compared them with Bagger-Gunion expressions ([3]) and §4§ one can fix errors made in Ref.[l|. 
Namely, ignoring error/misprint (instead of ir should be it 2 ) we see that the function D(s,u) does not coincide with 
A(s,tt) from Eq.Q: dependence on the subprocess invariants s and u are wrong, and normalization is also incorrect, 
because for the symmetric distribution amplitudes 
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In other words, the results of Ref. [l| contain additional numerical factor 1/4. 

In Sec. Ill authors tried to compute the subprocess cross section in the RC method. Explaining the choice of the 
renormalization and factorization scales made in their work, they write: "As is seen from (2.11) [Eq.( I14l) ]. in general, 
one has to take into account not only dependence of a{Q\ 2 ) on the scale Q\ 2 , but also an evolution of $(x, Q 2 ^) with 
Q\ 2 - The meson wave function evolves in accordance with a Bctc-Salpctcr-type equation. Therefore, it is worth noting 
that, the renormalization scale (argument of a s ) should be equal to Q\ — ~xu, Q\ — sx, whereas the factorization 
scale [Q 2 in $(x, Q 2 ] is taken independent from x, we take Q 2 = ■ In the frozen coupling approximation, as is 
evident also from the Bagger-Gunion expression (U), the choice for the factorization scales is s/2 and —u/2, whereas 
in the RC method without any reasons authors put n 2 F1 = ii 2 F2 = "Pt ■> Pt being the meson M transverse momentum 
square. Of course, such choice is not legitimate in the context of RC method. Actually authors did not compute 
the subprocess cross section using the RC method, but "generalized" their wrong expression obtained in the frozen 
coupling approximation by this way [the first line in Eq.(3.5) of Ref.[l|]: 
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where A = 1 — x (as it follows later from Eq. (3.7)). 

In our work [2l| it was proved that within the RC method the cross section of the subprocess jqi —> Mq2 is given 
by Eqs.( [6]), (JT)). Simple generalization of the Bagger-Gunion formula Q for the pseudoscalar and longitudinally 
polarized vector meson in order to use it in the RC method leads to wrong results, because in the RC method even 
for mesons with symmetric distribution amplitudes relations (1101) do not hold. The latter conclusion was made in Ref. 
[2lj and proved by explicit calculations, results of which were presented in Eqs. (30) and (32) of this paper. 

As is seen from Eq. (fl"6]) , in Ref. [l[ even for the renormalization scales (Ijq, /i^ 2 the choice [i 2 Rl = n 2 R2 = (1 — x)-p\ 
was made, which contradicts to principles of the RC method and declarations of the authors. Therefore, Eqs. (3.9)- 
(3.13) are wrong themselves and do not describe situation with the two scales. The Borel transforms presented in 
Eqs.(3.14)-(3.17) and resummed expressions shown in Eqs.(3.18)-(3.21) are also wrong in the case of two scales. We 
think that we can stop at this point our analysis of the "investigation" performed in Ref.[l] . 

Unfortunately, exactly the same errors were made by these authors in their paper published in Phys. Rev. D80, 
016003 (2009) and entitled "Infrared renormalons and single meson production in proton-proton collisions" (22j. 
Despite declarations of the authors the process that they tried to considere is not " single meson production" 

p + p -)• M + X, 

but "prompt photon and meson production" 

p + p^ 7 + M + X. 

In Eq.(2.10) of Ref. 22] authors present the diffential cross section of the HT subprocess qiq 2 — > 7t + (tt~) + 7 
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Here Q\ = — (1 — xi)u and Q 2 = —x\t. 

In Sec. Ill of Ref. [22| authors tried to compute the subprocess qiq 2 — > 7r + (7r~) + 7 cross section in the context of 
the RC method. But instead of calculation of the Feynman diagrams of the subprocess within the RC method they 
"generalize" the expression (1171) obtained in the frozen coupling approach. The correct expression for the subprocess 
cross section in the RC approach can be obtained from our formula ([6]) for the pseudoscalar mesons, after exchange 
of the Mandelstam invariants and by taking into account numerical factors. 

In this section the authors write:" As is seen from (2.11) [Eq.(fT8"|)], in general, one has to take into account not only 
dependence of a(Q\ 2 ) on the scale Qf 2 , but also an evolution of $(x, Q\ 2 ) wrtn Q\ 2- The meson wave function 
evolves in accordance with a Bete-Salpeter-type equation. Therefore, it is worth noting that, the renormalization scale 
(argument of a s ) should be equal to Q\ = (xi — l)u, Q\ — —txi, whereas the factorization scale [Q 2 in $(x, Q 2 } is 
taken independently from x; we take Q 2 — p\ " ■ Such choice for the factorization scales, i.e. equating them to the 
pion's transverse momentum square p^, in the HT subprocess calculations is wrong. The factorization scales had to 
be chosen equal to — u/2 and — t/2 in accordance with the standard prescriptions. As is seen from Eq.(3.5) of [22[ 
and the text above of the Eq.(3.7), not only the factorizations scales, but also the renormalization scales were chosen 
in a way that contradicts prescriptions of the RC method and statements of the authors. Indeed, in the first line of 
Eq.(3.5) having corrected the misprint (should be x instead of x\) we find 

Dm = ei f T ^W)*"(*.Q 2 ) + e2 u f 1 dx ^ Q y MM 2 )^ 



jo 1 x J X — j, 

where A = 1 — x. It is evident that in this case n R1 = [i 2 R2 = (1 — x)p 2 n 1 and they are not equal to Q\ = (xi — l)u, Q 2 — 
—tx\ as write authors in the paper. In other words, in Ref. [22j the factorization and renormalization scales were 
expressed using the transverse momentum square of the final particle, but not the momentum transfer carried by 
hard gluon. All computations in Ref. [22| were carried out with these incorrect choices for both the factorization and 
renormalization scales, which led to wrong expressions. 

Ignoring for the moment that expression for the function D(Q 2 ) p^]) in the RC method is wrong and does not 
describe pion production, let us nevertheless explane the correct approach to computation of such integrals. After 
correct choices for the scales it should have the form: 

D(i, u) = ei t f dx -s[-(l-^M^~u/2) + e2{i f dx a s [-xt\* M (x,-t/2), (20) 
Jo 1 — x Jq 1 — x 

Having expressed the running couplings a s [— (1 — x)u] and a s [— xi\ in terms of ct s (— u), a s (—t) with the aid of the 
renormalization group equation we find 

a{XQ ] = T+laX/t 

where a = a s /ir and t = An / f3oa s (Q 2 ) (there is an error/misprint in corresponding Eq.(3.4) of Ref.|22j|). Having used 
this expression the integrals in Eq. (|20[) can be recasted into the following forms: 
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and functions A[u, a s (—u/2)] , B[u, a s {— 1/2)] contain IR renormalon poles, number and locations of which depend on 
the meson (pion) distribution amplitude. These functions depend also on the factorization scales through a s (—u/2) 
and a s (—t/2), respectively. In other words correct treatment leads to expressions with two scales and two integrals, 
but not to formulas like (3.10)-(3.13) in Ref.[H| 
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t 4?r 
/3oa s (p|)' 

Therefore Eqs.(3.10)-(3.13) in Ref. |22j are wrong. Authors state that they found the relevant Borel transforms 
(3.10)-(3.13) of the resummed expressions . In the case of correct choice for the scales Eqs.(3.14)-(3.17) are wrong, 
because as it follows from our analysis and from Eq. (|21[) the corresponding Borel transforms have two components: 

B 1 [D(t,u)](u) = A[u,a s (-u/2)}, B 2 [D{t,u)](u) = B[u,a s (-t/2)]. (22) 

As a result, the resummed expressions should have also two pieces. Then, for example, for Eq. (l2T1) we get 

[D(t,u)r s = AV *l fneit A[Li(\V/\?, a.i-u/2)] + B[Li{\^)/\^ a s (-t/2)], (23) 

Po Po 



where 
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and the logarithmic integral Li (A) for A > 1 is defined in its principal value 

f x dx 

Li(X) = P.V. 
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(in Eq.(3.22) of Ref.[22| there is evident error/misprint). Consequently, Eqs.(3.18)-(3.21) presented in Ref.[22| are 
wrong as well. Our analysis of the Borel transforms and the resummed expressions presented in Eqs. (|20p - (|2"3"]l arc 
valid also for the first process 7 + 7— >Af + X [l|. 

We think that after this analysis the remaining part of the works [l[ and [13] , namely the leading twist contribution 
to the process and numerical calculations do not deserve further discussions. Unfortunately, all these fatal errors, 
wrong statements, numerous misprints were overlooked by referees of such prestigious journal like Phys. Rev. D. 
Therefore, we decided to inform the HEP community on "methods" used by these authors in their "investigations". 
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